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We perform a systematic study of the possible molecular states composed of a pair of heavy 
mesons such as DD, D* D, D* D* in the framework of the meson exchange model. The exchanged 
mesons include the pseudoscalar, scalar and vector mesons. Through our investigation, we find that 
(1) the structure X(3764) is not a molecular state; (2) There exists strong attraction in the range 
QQ , r < 1 fm for the D* D* system with J = 0, 1. If future experiments confirm Z+(4051) as a loosely 

■ bound molecular state, its quantum number is probably — 0^. Its partner state may be 

' searched for in the vr^'xci channel; (3) The vector meson exchange provides strong attraction in the 

CNj , D* D channel together with the pion exchange. A bound state solution exists with a reasonable 

cutoff parameter A ~ 1.4 GeV. X(3872) may be accommodated as a molecular state dynamically 
although drawing a very definite conclusion needs further investigation; (4) The B* B molecular 
state exists. 
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PACS numbers: 12.39.Pn, 12.40.Yx, 13. 75. Lb 



Q_^• I. INTRODUCTION 

' Since the observation of Z+(4430) [l|, the Belle Collaboration reported two new resonance-like structures ^+(4051) 
. and ^+(4248) in the tt+XcI niass distribution in the exclusive K^ir^Xci decay. Their masses and widths are 

mz+(405i) = (4051 ± 141^5) MeV, Tz+(m5i) = (821^1+47) MeV and mz+(4248) = (4248± 14t^4+i80) ^eV, r^+(4248) - 
(ITTlgglg^®) MeV 0. These charged hidden charm signals are good candidates of either tetraquark states or heavy 
^ molecular states if they are confirmed by future experiments. 

, Recently the EES Collaboration reported an anomalous line-shape observed in the range of 3.650 GeV to 3.872 
' GeV by analyzing the cross section of e^e^ — s- hadrons. The anomalous line-shape is composed of two possible 
\ enhancement structures around 3.764 GeV and 3.779 GeV respectively [3]. The later one is consistent with the well- 
• ' established ■0(3770) while the mechanism of the first structure is not clear now, which is denoted as AT (3764) in this 
\ work. 

In the past five years, a series of observations of the charmonium-like AT, Y, Z states, especially those states near 
the threshold of two charmed mesons have stimulated the interest in the possible existence of heavy molecular states 
greatly. The presence of the heavy quarks lowers the kinetic energy while the interaction between two light quarks 
could still provide strong enough attraction. 

In fact, Voloshin and Okun studied the interaction between a pair of charmed mesons and proposed the possibilities 
5^ ' of the molecular states involving charmed quarks more than thirty years ago ,4j. De Rujula, Geogi and Glashow 
\ speculated ?/'(4040) as a D*D* molecular state [5j]. Tornqvist studied the possible deuteron-like two-meson bound 
states such as DD* and D*D* using the quark-pion interaction modelteL Dubynskiy and Voloshin suggested the 
possibility of the existence of a new resonance at the D*D* threshold Zhang, Chiang, Shen and Zou studied 

the possible S-wave bound states of two pseudoscalar mesons with the vector meson exchange in Ref. [9]. The 
experimental observation of .A (3872) and Z+(4430) motivated the extensive discussion of X( 38 72) as a DD* molecular 
state [l3,[il|,[il,[il,[il,[il,[il[i3,lll and Z+(4430) as a D[D*{DiD*) molecular state [ll, H, Ell, [H . 

It's interesting to note that the ^'"'"(4051) enhancement announced by the Belle collaboration and the X(3764) 
signal reported by the BES collaboration are near the thresholds of D*D* and DD respectively. One may wonder 
whether they could also be candidates of heavy molecular states. In this work we perform a systematic study of 
three types of possible heavy molecular states: the DD/BB system (P-P), the D*D* /B*B* system (V-V), and the 
D*D/B*B system (P-V) using the formahsm developed in Ref. [Tsl. [2ll. [2a| . 
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This paper is organized as follows. After the introduction, we introduce some notations and present the flavor wave 
functions of the S-wave P-P, V-V and P-V systems constructed by two heavy flavor mesons. In Sec. Illli we collect 
the effective Lagrangians in the derivation of the effective potential. Sections IIVIIVII are for the P-P, V-V, P-V cases 
respectively. The last section is the discussion and conclusion. 

II. FLAVOR WAVE FUNCTIONS OF HEAVY MOLECULAR STATES 

We study the possible molecular states composed of two pseudoscalar (P-P) heavy mesons, two vector (V-V) heavy 
mesons, one pseudoscalar and one vector (P-V) heavy mesons. The masses of the — 0^, 1^ heavy mesons are 
taken from PDG [H and collected in Table HI 







charmed-up(down) meson 


charmed-strange meson 




charged 


mass (MeV) neutral mass (MeV) 


charged mass (MeV) 


0" 

1" 




1869.3 D" 1864.5 
2010.0 D*° 2006.7 


Df 1968.2 
D** 2112.0 






bottom-up(down) meson 


bottom-strange meson 


r 


charged 


mass (MeV) neutral mass (MoV) 


neutral mass (MeV) 


1" 




5279.1 B" 5279.5 
5325.1 B*° 5325.1 


5365.1 
Bf 5412.0 



TABLE L The masses of heavy mesons in the H doublet [23l |. 



The P-P type is categorized as two systems, i.e. V — V, B — B. Here V, V, B and B denote , , Df), 
{D°,D- ,D-), {B+,B^,B°) and {B-,B°,B^) triplets respectively. In the following, we illustrate their flavor 
wave functions with the V — V type as an example. Since charmed mesons belong to the fundamental represen- 
tation of flavor SU{3), the V — V system form an octet and a singlet: 3 ® 3 = 8 © 1. The corresponding fla- 
vor wave functions are = D^D+, $+ = D°D+, $o = D-D+, $" = -1={D°D° - D-D+), l-^ = D-D+, 

$8 = ^(D^D° + D-D+ - 2D-D+), <p- = D-D", <i>- = D-D^, $i = -j={D^D" + D- D+ + D-D+). Through 

ideal mixing, we have = ■^(D'^D" + D^ D^) and ^^j^ = D~Df, where $g is an isoscalar while $° is an isovector. 

The above states are shown in Fig. [T]and collected in Table HTl together with the B — B type states. The flavor wave 
functions of the D*D* and B*B* systems are constructed similarly, which are denoted as $** and respectively. 



$0 $+ 




FIG. 1: The molecular multiplets composed of charmed and anit-charmed mesons. 



We label the T> — T>* and B — B* systems as $* and ft* respectively and list their flavor wave function in Table 
Hm Here we need distinguish the C parity for <^*°, and $*? in the V-V* system and and fl*^ in the 

B — B* system. The parameter c — ^1 corresponds to C = ±1 respectively as pointed out in Ref. [3, [2l|. In this 
work we label those states with the negative charge parity with a hat such as <I>g°, $*°, fi***, jig" and 51*°. The 
charge parity of X(3872) is positive. 



3 



V-V B-B 



State 


wave function 


State 


wave function 


$+ 

^ S 

$- 
$7 

3>?i 


D°D+ 

D-Dt 
-^{D D — D D^) 

D-D+ 

D-D" 

D-D'' 
^(D°D° + D-D+) 

D-Dt 


n+ 

nO 

Q- 


B°B° 
-^[B^B —BB) 



TABLE II: Tlie flavor wave functions of the V — D and B — B systems. 

















State 


wave function 


State 


wave function 


$:+/$:+ 




^{D-'^Dt + cD''Dl+) 










$*+ /$*+ 




-^{D'-°D+ +cD°D''+) 














^{D*-Dt +cD-D:+) 


nf/nf 












- D*-D+) + c (_D''_D*'' - _D-D*+)] 






B- - +c(B+B*- 


-B^B'")] 












^(BfB«+cB°B*°) 




$— /$*- 






Q*-/Ci*- 




■j={B*°B- +cB°B*-] 




*r/-i'r 










j^{BfB-+cB°B*-] 








+ D—D+) + c {D°D*° + D-D*+)] 


nf/nf 




-B- +B^''B°)+c{B+B*- 


+ B°B*°)] 






^{DrDt +cD:D*+) 






^^{B'lB'l+cB'lBf) 





TABLE III: The flavor wave functions of the T> — D* and B — B* systems. 



III. THE EFFECTIVE LAGRANGIAN IN THE DERIVATION OF THE POTENTIAL 



In the derivation of the potential, we need the the effective Lagrangians, which are constructed based on the chiral 
symmetry and heavy quark symmetry: psi . [26l . [27j 



where the field H is defined in terms of the (0,1) doublet 



A'^i^ is the axial vector field with definition 



with ^ ~ exp(zP//7r) and = 132 MeV. The octet pseudoscalar and nonet vector meson matrices are 



(1) 



(2) 



(3) 



/ Jll + JL 

77- 



(4) 



4 



V 



V K*- K*° 




The effective interaction Lagrangians at the tree level from Eq. ([T]) are 

CvVcr = -2 7711) gaVaVla, 

r^.^.v = zgw'V {Vl'^^d^Vl^b - Vlbd^Vl''^) {N^)ab + ^ifv'V'yDl\Dlb{d''N'^ - d''W)ab, 
Cv^va = 2vav^g,VTV*^^o, 
C-n.vr = -igv'Vv{DaV*X-Vl,vl)d^'Tab. 

where I?(*) = ((L)0)(*), {D^)^*\ (D^ )(*'). The relevant coupling constants are 

gv'VT 2g (3gv 
gvvv — , = — -7-, gvvY — gvv 



y/mvmv /tt V2 

f fvVY ^gv _ ™p _ 5Tr 

Jx>*r>v — —7=-, gv — 9<y - 7-7=, yo) 

mv V2 U 2V6 

where , P and A are parameters in the effective chiral Lagrangian that describe the interaction of heavy mesons 
with hght vector mesons Following Rcf. [H, we take g = 0.59^^= 0.9 and A = 0.56 GeV"^ g^ = 3.73 [H]. 

In this work, we adopt the same formalism developed in Refs. [l^, [U, [l^] to derive the effective potential. We 
compute the amplitudes of the elastic scattering of two heavy mesons using the above effective Lagrangians. In order 
to account for the structure effect of every interaction vertex, we introduce the monopole type form factor (FF) [1, [131 

Here A is the phenomenological parameter around 1 GeV, and q denotes the four-momentum of the exchanged meson. 
The FF also plays the role of regularizing the potential by imposing a short-distance cutoff to cure the singularity of 
the effective potential. Then we impose the constraint that the initial states and final states should have the same 
angular momentum. After averaging the potentials obtained with the Breit approximation in the momentum space, 
we finally perform Fourier transformation to derive the potentials in the coordinate space. In the following sections, 
we illustrate the effective potentials in detail for the P-P, V-V and P-V systems. 

IV. THE V-V CASE 

A. The potential of the P-P system 

The quantum number of the S-wave T) — T) system is — 0+ and the C parity of the neutral states is positive. 
In this work we consider the possible molecular states bound by the force from the light pseudoscalar meson, vector 
meson and scalar exchange. Parity and angular momentum conservation forbid the exchange of a pseudoscalar meson 
between the V — f) pair. No suitable meson exchange is allowed for the and $°($^) cases. For the other states, 
the effective potentials in the momentum space from the vector meson and a meson exchange are 



1 , q' 



2 
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where and my denote the masses of charmed meson and exchanged vector meson respectively. /^^I'^i and /J^^f 2 
arise from the coefficients related to the exchanged meson, which can be read from the nonet vector meson matrix in 
Eq. ([5]). After making the Fourier transformation, the three independent structures in Eqs. ([8])-(l9]) read as 



where 



Y[A, m, r] 
Z[A, m, r] 
X[A,TO,r] 



q2 ^ ^2 



A(A,m,r), 
>y(A,m,r), 



1 



1 d 



-Ar 



SttA 



or V or ' 



1 9 



d 



Z[A, m, r] 



(10) 
(11) 

(12) 



(13) 
(14) 

(15) 



with ^ = VA^ — m^. We have adopted the monopole FF in Eq. ([7]) to regularize the potential in Eqs. ([8])-((9]). Now 
the effective potentials in the coordinate space read 



Vv(0 



tYViVi JYV2V2 
'1 J -'2 



^9 



y[A,mv,r] + 



A [A, my, r] 



V,(r) = -g2y[A,m,,^ 
The exchange potential for $^ in coordinate space is 

1 1 



(16) 
(17) 



nr)Totai - Vv(r) 



V2' ^/2 



Vv(r) 



1 1 



+ V.(r) 



2 



Y[A,m^,r] Y[A,mp,r] A[A,TO„,r] , A[A,TO^^,r] 



- glY[A,m^,r]. 



2 2 

The symbol means the SU{2) symmetry is assumed. For the and <i>g states, the exchange potentials are 

1 1 



(18) 



V(.)?rP = i Vv(r) 



-Vv(r) 



+Vv(r) 



9vw I 



1 1 

1 1 



=F 2Vv(r) [1, l,TO£,o,m£,+ ,mp±] 
1 1 



Vv(r) 



V2' V2 



(1t2) y [A, mp,r] + 



A[A,TOp,r] 



Y[A,m^,r] + 



X[A,', 



^g2y[A,m^,r]. 



(19) 



Clearly the potential of $0 is the same as that of due to the SU(2) symmetry. The effective potential of <i>°]^ is 



^i'^)Totai = Vv(f) [1, l,mD,,mD,,m^ 



-9vvy >"[A,TO0,r] + 



A[A,m0,r] 
4m|j rn^ 



(20) 



where only the meson exchange is allowed. 
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B. Numerical results for the P-P system 



The input parameters include the coupHng constants in Eq. ^ , the masses of heavy mesons hsted in Table HI and 
the exchanged meson masses = 139.5 MeV, nip — 775.5 MeV, = 600 MeV, m,, — 547.5 MeV, = 782.7 
MeV and = 1019.5 MeV. 

We first plot the variation of the effective potentials of the T) — T) system with the cutoff A = 1 GeV in Fig. [2l For 
the <i>='=(<i>°) state, there exist the p, w and a meson exchanges. The p exchange force is repulsive while both the lo 
and a meson exchange forces are attractive. The p and lo exchange forces cancel each other almost exactly. Thus the 
total potential is attractive. For the state <I>gj, the </> meson exchange potential is attractive. The exchange potentials 
of p, LO and a mesons all are attractive for <I>g. Thus its total effective potential is attractive. 




FIG. 2: Diagrams (a), (b) and (c) give respectively the variation of potentials of $ ($)", and $g states on r. 



We adopt the MATSLISE package to solve the Schrodinger equation with the effective potentials of <i>='=(<i>°), <i>g]^ 
and <i>g states. MATSLISE is a graphical Matlab software package for the numerical study of regular Sturm-Liouville 
problems, one-dimensional Schrodinger equations and radial Schrodinger equations with a distorted coulomb potential. 
It allows the fast and accurate computation of the eigenvalues and the visualization of the corresponding eigenfunctions 

The binding energies with typical values of A for $^($°), $2]^ and $g are presented in Table HVl if there exists the 
solution by solving Schrodinger equation, rmis denotes the root-mean-square radius with the unit of fm. For 
we can not find bound state solutions in the range of A < 10 GeV, which indicates <I>='=(<I>°) does not exist. For both 
<^1i and $g, there exist bound state solutions with A around 0.5 GeV. When the binding energy E becomes smaller, 
the r^rns becomes larger. Similar observations hold for the B — B system. Now the reduced masses of VP^i and fig are 
heavier. Hence the kinetic term is smaller. 51"]^ and Jig states appear with a larger A around 1.1 GeV. As an example, 
we show the dependence of the binding energy of <I>°j^ on A in Fig. [31 Moreover, we increase and reduce both the 
vector coupling constant and scalar coupling constant by a factor of two to see the variation of the binding energy 
and the cutoff parameter. The munerical results are collected in Tables fVl and IVll 



V. THE V* - V CASE 



A. The potential of the V-V system 



The possible quantum numbers of the S-wave V-V system are = 1+, 1+, 2+. The C parity is + for the neutral 
states. The exchanged mesons include the pseudoscalar, vector and a mesons. The exchange potential reads 



Top 



1 2 
^ yvD 



(21) 
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V-V 


States 


A (GcV) 


E (UvN) 


ri-m„ (fm) 


$± (<!,") 






0.58 


-10.14 


1.84 




0.60 


-3.46 


2.80 




0.50 


-23.42 


1.44 




0.53 


-4.93 


2.52 




B - B 


State 


A (GeV) 


E (MeV) 


rrms (fm) 












0.70 


-9.24 


1.84 




0.72 


-1.78 


2.86 




1.10 


-7.13 


1.95 




1.15 


-23.34 


1.23 



TABLE IV: The bound state solutions for and <il and fi^i and fil 




0.55 0.56 0.57 0.58 0.59 0.60 



A (GeV) 



FIG. 3: The dependence of £ on A for 



for the pseudoscalar meson exchange and 



V|/l(q) 



"-'l -'2 



1 



C[J] ^ + gl,^,^C[J] ^— ^ 



Tod* TOD* 



(22) 







State 


A (GeV) 


E (MeV) 


rrms (fm) 












0.72 


-23.00 


1.30 




0.74 


-8.16 


1.90 


<J>^ 


0.60 


-37.20 


1.17 




0.62 


-11.39 


1.77 




B-B 


State 


A (GeV) 


E (MeV) 


rrms (fm) 


n*(n") 










0.82 


-16.70 


1.44 




0.83 


-7.69 


1.92 




0.91 


-11.84 


1.66 




0.92 


-21.27 


1.32 



TABLE V: The bound state solutions for ^^(^O), $2i and and Q.^{QP), figi and nl if we increase the coupling constants 
by a factor of two. 
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State 


A (GeV) 


E (MeV) 


rrms (fm) 


$±((1)") 








TTi 


0.38 


-24.49 


1.50 




0.39 


-18.89 


1.63 


$8 


0.38 


-15.38 


1.78 




0.39 


-10.61 


2.02 




B-B 


State 


A (GeV) 


E (MeV) 


rrins (fm) 


V / 










0.50 


-30.69 


1.29 




0.52 


-17.64 


1.55 




0.46 


-29.24 


1.37 




0.48 


-14.15 


1.73 



TABLE VI: The bound state solutions for $*($°), t&^i and $8 and $1^(51"), and S7g if we reduce the coupling constants 
by a factor of two. 



for the vector meson exchange. For the cr exchange, 



vW(q) = -glC[J] 



q2 + m?. 



with 



A[J\ 
B[J] 
C[J] 



(lAi;lA2|J,m)(lA3;lA4|J,m)^ [e^ • {qxel^^^' ' (q x ef^) 



Ai A2 A3 A4 



(lAi; IA2I J,m)(lA3; IA4I J, m)^^ [(ei ' q)(e2 • ^M ' e^) + (ci-s)], 

Ai A2 A3 A4 

Y (lAi;lA2|J,m)(lA3;lA4|J,m)(et^ -e*^^) 



A2 , *A,_ 



(23) 

(24) 
(25) 
(26) 



Ai A2A3A4 

where ei is the polarization vector, and c.t.s denotes the cross-terms, i.e. — (ei • q)(e4 • q)(e2 ■ £3), (£3 ■ l){^4 ■ q)(ei ■ ^2) 
and —(£3 • q)(e2 • q)(ei ■ £4)- When considering the different systems with — 0"^, 1+, 2+, we impose the constraint 
on the scattering amphtudes that the initial states and final states should have the same angular momentum. Then 
we average the potential in the momentum space, i.e., making the substitutions q^yz^ 9^/3. The coefficients of 
A{J), B{J) and C(J) are listed into Table |VlIl 



J 


A{J) 


B{J) 


C(J) 





2/3 


4/3 


1 


1 


1/3 


2/3 


1 


2 


-1/3 


-2/3 


1 



TABLE VII: The values of A{J), B{J) and C(J) for the cases of = 0+, 1+, 2+ systems. 



After Fourier transformation, we get the potentials in the coordinate space 



— ^1 ^2 



1 



gl'v-rh ' ' h ^-4[J]Z[A,mp,r], 



Tzf^ C[J] A[A, mv, r] + gl,T,,^C[J] F [A, my, r] 



m-D'Tn-D* 



■ B[J] Z[K, my, r] 



(27) 



(28) 



V.(r)W = -glC[J]Y[K,m,,r]. 
The total effective potentials for the and states are 



-^,--^,mD'0,mD',mn 



-.,mD,+ ,mD^,m, 



where only the 77 meson exchange is allowed. Considering the SU(2) symmetry, one further gets 



For the state, the potential is 



1 



+ "14/' (r) [1,1, rriDj , mcj , m,^] 



C[ J] X[A, m^, r] + ^l^.p-v C[J] 1^[A, m^, r] 



+ 



■S[J]Z[A,m^,r] 



For the exchange potential reads 



'^V) Total 



1 1 

71' 71 
1 1 

7i'~7I 



m£).o,m£).+ ,m,r 



+ V[/I(r) 



1 1 

—^,mn'0,mn'+,mn 

1 1 



72' 72 



, m£).o, m£).+ , ma 



Z[A,m^,r] _ Z[A,m^,r] 
8 24 



1 

+ 2 



-5^C[J]y[A,m,,r] 

1-2 



x[A, ^ 5l,.2,.v C[J] Y[A, m„ r] + g^jj ^j^, m„ r] 



9V'V'Y 



C[J] X[A, m^, r] + gl,-!,,^C[J] Y[A, m^, r] + 1/%^ B[J] Z[A, m„, r 
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The potential for the state is 



JTotal 



[J] 



1 1 



\/2' V2 
1 1 

1 1 

1 1 



if2Vp\r) [l,l,m£,.o,TO_D.+ ,TO^±] 
1 1 



V6' V6 



1 



1 



+ V|/l(r) 



1 1 



1 1 



^/2'^/2 
(IT 2) 



2VW(r) 



Z[A,TO^,r] 
24 



(IT 2) 



9v*T>*Y 



4/: 



C[J] X[A, mp, r] + gW^^CiJ] Y[A, m^, r] + -^1^ B[J] Z[A, nip, r 



C [J] X [A, m„ , r] + .g|,.p. V C [J] F [A, , r] 



4/ 



^e[J]Z[A,m„,r] 



(35) 



where =F corresponds to $**"($**"') respectively. As a cross-check, the potential for $**"' is the same as that for 
from the SU(2) symmetry. For the V-V systems, we need only consider four independent types of potentials 
corresponding to [<^T^, ^**° , ^**% [$**±, and listed in Eqs. respectively when 

we consider the SU(2) symmetry. For states with different quantum numbers, we simply change the values of A[J]^ 
B[J] and C[J] in Table EIll 



B. Numerical results for the V-V system 



The quantum numbers of the S-wave V-V system are = 0+,l+,2+. The variation of the oscillating effective 
potential of $**±($**o^ ^**0) ^-^^-^ 

r is shown in Fig. [H where only the rj meson exchange is allowed. The effective 
potentials of and $1*° with = 0+, l"*", 2+ are presented in Fig. where the exchanged mesons include 

TT, ry, p, Lo and a. For its effective potential is shown in Fig. [71 where both the i] and meson exchange is 

allowed. 
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FIG. 4: The effective potential of $***(*r°, *r°)- 
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FIG. 5: The effective potentials of $**='=($**'') with different quantum numbers. The thick solid line is the total effective 
potential. 





With the above effective potentials, we solve the Schrodinger equation to find the bound state solutions for the V-V 
system. Numerical results are collected in Table IVini For the T)* — V* states with = 0+, 1+, a bound state exists 
only for A much larger than 1 GeV. Especially for ($**°, and <j)**±(-(j)**o^ ^^^^ jP ^ 2.+, the existence of 

a bound state requires the values of A be larger than 8 GeV. For $**±(<I>**o, and $**±($**0) with = 2+, no 

bound states exists for A < 10 GeV. The cutoff parameter A is a typical hadronic scale, which is generally expected 
to be around 1 '--^ 2 GeV. If A is much larger than 2 GeV or much smaller than 1 GeV in order to form a bound 
state, we tend to conclude that there does not exist a heavy molecular state for these systems. We find bound state 
solutions for$|*o and with = 0+, 1+, 2+. 

In Table IVIIIi we also give the numerical results for the B* — B* system. There do not exist molecular states 
and with = 2+ since we can not find bound state sojutions for A = - 10 GeV. 

Different from the V* - V* system, there exist possible molecular states with = 0+,l+ with 

more reasonable A close to 1 GeV. We also find the bound state solutions for S7g*° and ri*j['°. The reason is simple. 
The larger the reduced mass for B* — B*, the lower the kinetic energy, the easier to form a molecular state. With 
$** as an example, we illustrate the dependence of on A in Fig. [H] When all coupling constants are enlarged by a 
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factor of two, the numerical results are collected in Table HXl When all coupling constants are reduced by a factor of 
two, the results are collected in Table 1X1 




-120 I 1 1 1 1 1 1 1 1 

48 4.9 5,0 5 1 5,2 
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FIG. 8: The dependence of binding energy i5 on A for $3*. 











V* 
















J^ = 0+ 






J^ = l+ 






J^ = 2+ 




State 


A (GeV) 


E (MeV) 


Trms (fm) 


A (GeV) 


E (MeV) 


rrms (fm) 


A (GeV) 


E (MeV) 


r-rms (fm) 




4.80 


-5.26 


1.92 


8.70 


-1.60 


2.56 










5.00 


-39.42 


0.72 


8.90 


-21.21 


0.96 










4.00 


-5.21 


1.99 


9.80 


-8.40 


1.61 










4.20 


-19.02 


1.08 


10.00 


-16.53 


1.19 










0.52 


-18.32 


1.48 


0.50 


-37.61 


1.18 


0.52 


-11.71 


1.85 




0.53 


-8.61 


1.95 


0.52 


-15.03 


1.63 


0.50 


-24.01 


1.47 




0.63 


-12.91 


1.55 


0.60 


-17.27 


1.45 


0.54 


-23.64 


1.46 




0.64 


-7.13 


1.93 


0.62 


-6.56 


2.11 


0.56 


-13.01 


1.77 












- S* 














J^ = 0+ 






J^ = l+ 






J^ = 2+ 




State 


A (GeV) 


E (MeV) 


rrms (fm) 


A (GeV) 


E (MeV) 


rrms (fm) 


A (GeV) 


E (MeV) 


rrms (fm) 




2.30 


-0.57 


4.76 


3.80 


-3.98 


2.21 










2.40 


-14.90 


1.17 


3.90 


-17.91 


1.06 










1.50 


-10.17 


1.53 


2.70 


-3.60 


2.41 










1.60 


-21.75 


1.10 


2.90 


-15.55 


1.22 










0.57 


-23.16 


1.36 


0.57 


-32.70 


1.23 


0.61 


-18.07 


1.49 




0.58 


-6.42 


2.28 


0.59 


-5.65 


2.39 


0.63 


-8.39 


1.93 




0.74 


-12.12 


1.56 


0.70 


-33.5 


1.09 


0.64 


-35.37 


1.25 




0.75 


-4.48 


2.37 


0.72 


-12.68 


1.58 


0.65 


-27.12 


1.35 



TABLE VIII: The numerical resuhs for the V - V and B* - B* systems. 



VI. THE V - V CASE 
A. The potential of the P-V system 



Fig. [S] shows the the general scattering channels in the derivation of the exchange potential. Different from the 
P-P and V-V systems, we need consider two additional crossed scattering diagrams for the P-V system, i.e. Fig. [9] 
(b) and (c). 
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V 


-V* 














J- = 0+ 






J- = l+ 






J^ = 2+ 




State 


A (GeV) 


E (MeV) 


rrms (fm) 


A (GeV) 


E (MeV) 


rrms (fm) 


A (GeV) 


E (MeV) 


rrms (fm) 




1.90 


-16.07 


1.14 


2.90 


-14.48 


1.18 


- 


- 


- 




1.85 


-4.7 


2.05 


2.80 


-1.14 


3.74 


- 


- 


- 




1.00 


-5.42 


2.06 


1.80 


-6.28 


1.89 


- 


- 


- 




1.20 


-36.75 


0.92 


2.00 


-24.25 


1.04 


- 


- 


- 




0.58 


-45.86 


1.05 


0.60 


-15.53 


1.62 


0.65 


-32.94 


1.15 




0.59 


-13.02 


1.73 


0.61 


-3.64 


2.64 


0.66 


-27.10 


1.21 




0.76 


-39.20 


0.97 


0.74 


-38.98 


1.01 


0.70 


-23.96 


1.38 




0.78 


-7.90 


1.82 


0.76 


-12.73 


1.55 


0.72 


-11.31 


1.78 












- e* 














J^ = 0+ 






J^ = l+ 






J^ = 2+ 




State 


A (GeV) 


E (MeV) 


rrms (fm) 


A (GeV) 


E (MeV) 


rrms (fm) 


A (GeV) 


E (MeV) 


rrms (fm) 




1.20 


-9.26 


1.50 


1.60 


-6.27 


1.79 


- 


- 


- 




1.30 


-79.29 


0.60 


1.70 


-44.38 


0.73 


- 


- 


- 




0.50 


-17.96 


1.41 


0.80 


-7.74 


1.80 


- 


- 


- 




0.48 


-16.36 


1.48 


0.90 


-25.86 


1.11 










0.60 


-81.06 


0.89 


0.62 


-46.66 


1.15 


0.70 


-245.16 


0.50 




0.61 


-10.73 


2.15 


0.63 


-11.49 


1.98 


0.71 


-352.21 


0.49 




0.83 


-25.32 


1.17 


0.82 


-34.93 


1.06 


0.80 


-26.93 


1.27 




0.84 


-5.08 


2.18 


0.83 


-15.01 


1.47 


0.82 


-11.66 


1.69 



TABLE IX: The numerical results for the V* — V* and B* — B* systems after increasing the coupling constants by a factor of 
two. 











V 


- V 














J^ = 0+ 






J^ = l+ 






J^ = 2+ 




State 


A (GeV) 


E (MeV) 


rrms (fm) 


A (GeV) 


E (MeV) 


rrms (fm) 


A (GeV) 


E (MeV) 


rrms (fm) 


$:**(<&:*", 3>r) 


















0.40 


-18.40 


1.57 


0.40 


-13.09 


1.82 


0.39 


-10.14 


2.10 




0.41 


-12.23 


1.81 


0.41 


-8.43 


2.12 


0.40 


-6.67 


2.42 




0.46 


-6.43 


2.21 


0.40 


-27.87 


1.37 


0.40 


-12.05 


1.95 




0.47 


-3.70 


2.76 


0.42 


-16.15 


1.63 


0.41 


-8.61 


2.18 










B* 


- B' 














J^ = 0+ 






J'' = 1+ 






J^ = 2+ 




State 


A (GeV) 


E (MeV) 


rrms (fm) 


A (GeV) 


E (MeV) 


rrms (fm) 


A (GeV) 


E (MeV) 


rrms (fm) 




6.80 


-9.33 


1.43 
















6.90 


-23.17 


0.93 














n"±(!^**") 


6.10 


-7.37 


1.64 
















6.20 


-14.64 


1.20 
















0.48 


-26.88 


1.31 


0.48 


-19.59 


1.51 


0.48 


-11.34 


1.92 




0.50 


-8.85 


1.97 


0.50 


-6.57 


2.25 


0.50 


-4.18 


2.75 




0.58 


-11.18 


1.68 


0.53 


-29.28 


1.23 


0.50 


-18.73 


1.60 




0.60 


-3.11 


2.85 


0.55 


-15.49 


1.54 


0.52 


-9.95 


1.97 



TABLE X: The numerical results for the "D* — T>* and B* — B* systems after reducing the coupling constants by a factor of 
two. 



In terms of the effective Lagrangian in Sec. IIIIl the a meson exchange potential in the momentum space is [3l| 
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; ^ 

V * V 
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FIG. 9: The diagrams in the derivation of the effective potential of the P-V system. 



The pseudoscalar meson exchange potential reads [l£ 
Vp(q) J2 , m-Di , mx). , mp 



1 



12tox>i™i3, 



yvvp-'-i ^: 



foi'TT 

for 77 



where = - ml, ^l' 



(37) 



Qq and go = — tod^. For the vector meson exchange, there exist two 



types of potentials from the direct and crossed scattering channels: 



gvvvgvv'v^i ^2 



1 



,(38) 



V^'-°^^(q) 



h J2 ,m-pi,mT,-,mY 



_ _8 2 j-YVi-D'Yv^vi q; 



q^ + C' 



(39) 



where Cv = y/m^^'--^. One notes that my is larger than qq for the P-V case. After Fourier transformation, we 
obtain the effective potentials in the coordinate space 



(40) 



f ^ \ j-PViv; j-pv;vi 



mp 



1 



^2 j-PViv; jPv;vi 



'Qv'DP-'-l 



U{r) forTT 
-Z[A,fi',r] for 77 



(41) 



-gvvY gv 



-VY J-i ^2 1 



, TOy, r 



+ r[A,TOv,r] , 



(42) 



(r) h ^2 ,mv„mv'^,mY 



Here 



U(r) 



47rr 



[cos{fir) 



» .2 j-vv.v; vv;v, , 



/32 a 



with /3 = - ml and a = v/A2 - q^. The total effective potentials of $*"/$*", $f/i>f 



*3 



For $*±/(i)*±^ the total effective potential reads 



g^'vp 



:Z[A,//,r] 



^('')Tota/ = l^gvVYgV'V' 



+ y[A,mp,i 



■^gvvYgv'v 



(43) 



(44) 



+ F[A,TO^,r] 



24m£iTOD* 



U{r) 



gV'VP /zfA ,,/ 1 4 2 ,7r. , „1 1 4 J-2 



(45) 
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For $*"/$*o, and $8°/$;°, we have 



1 



X[A, mp,r] 



+ r[A,mp,r] 



1 



X[A,m^,r] 



2 2 

gvvv U{r) — ^T^'T^^ 



72m]jm]j 



-Z[A, fi', r] - .^yZ[A, Cp, r] + ^/|,.pvZ[A, C^, r] | 



(46) 



V(r) 



Tota/ ~ 



-gvvYgvv 



Y[k,m^,r\ 



c-f^ 



IStoxJjTO^i. 



-Z[A,Ai',r], 



(47) 



)Total 



-gwygv'v 



AmDiTi]j*m'^ 



+ Y[A,mp,r] 



1 



-gwygw* 



X[K, 



UP 



5r"r>p 



4 



+ y[A, Too,, r] 



^g^r[A,m^,r] 



(48) 



The parameter c is =Fl for the P-V systems with the positive and negative charge parity respectively. As expected, 
the potential of ^*'^ and arc the same as that of and respectively. 



B. Numerical results for the P-V system 

We show the variation of the potential with r in Fig. [TOl The total effective potential for $*='=($*0) and <I>g° is 
attractive while the potential of $*=•=($*") and ^g*' is repulsive with A = 1 GeV. Especially the neutral states with 
positive charge parity are very interesting. The only candidate of AT (3872) with the attractive potential is $1". We 
collect the numerical results in Table IXll There do not exist molecular states $^($*'^), ri='=(r2*°) and (l^'^ because 
of the strong repulsive potential. Bound state solutions exist for ($*°, $*") and il*='=(r2*'', f]*"), l]*='=(f2*'^) only 
with a very large A. We show the dependence of the binding energy of $2i on A in Fig. [TT] Moreover we collect the 
numerical results in Table IXU when all coupling constants are increased by a factor of two. When all the coupling 
constants are reduced by a factor of two, the results are shown in Table IXlTTl 



VII. DISCUSSION AND CONCLUSION 



We have systematically studied the P-P, V-V and P-V systems, which is composed of a pair of heavy meson and 
anti- meson. We summarize our numerical results from Sees. IIVII VII in Tables |XIV| IXVI and IXVll We use the symbols 
'O' and '{^' to indicate whether a molecular state exists or not with the coupling constants in Section WUl In the 
following two cases, we tend to label the system with the symbol '(3)': (1) there does not exist a bound state at all 
because of the repulsive potential; (2) a bound state solution exists only with either A > 3 GeV or A < 0.5 GeV. If a 
bound state exists when 0.5 < A < 0.9 GeV or 2.0 < A < 3.0 GeV, we mark this state with '?'. In this case the result 
is so sensitive to the cutoff parameter that we can not draw a definite conclusion. Only when a loosely molecular state 
exists with 0.9 < A < 2.0 GeV, we tend to label the state with the symbol 'Q'- We want to emphasize the above 
criteria may be biased due to the authors' personal judgement. The readers are encouraged to consult the numerical 
results in Sections HVUVII and draw their own conclusions. 
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FIG. 10: The exchange potential of the P-V system with the typical value A = 1 GeV. The thick solid line is the total effective 
potential. 




It's interesting to compare our results with those near-threshold structures observed recently. 
1. Z+(4051) 
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of ill 


A (GeV) 


E (MeV) 


rrms (fm) 


<P (T> (T> 


- 


- 


- 


^*+ ^*() 
q>*^ Q*^ qp*^ 


8.90 


-11.74 


1.29 




9.00 


-25.38 


0.88 




10.5 


-12.86 


1.35 




- 


- 


- 


-.- -fci 1 

$ 1 


0.55 


-19.51 


1.55 




0.57 


-10.11 


1.93 


^ si 


0.60 


-17.96 


1.43 




0.62 


-6.94 


2.01 




0.52 


-13.02 


1.79 




0.54 


-5.25 


2.38 


^8 


0.50 


-28.61 


1.30 




0.52 


-9.42 


1.91 




B-B' 




A (GeV) 


E (MeV) 


rrms (fm) 


A*± o*'* rS*'* 








(~\ * zt (~\ * U (~) * U 


3.80 


-3.07 


2.20 




3.90 


-20.23 


1.00 




2.80 


-4.20 


1.41 




3.00 


-21.13 


0.96 










^iil 


2.10 


-5.94 


1.99 




2.20 


-9.79 


1.65 




0.72 


-12.40 


1.59 




0.73 


-5.95 


2.08 




0.92 


-10.19 


1.13 




0.94 


-18.67 


1.03 







TABLE XI: The possible bound state solutions for the P-V system. 



The charged broad structure Z+(4051) was observed by Belle m the tt+XcI. Its central mass value is slightly 
above the D*D* threshold. First we want to emphasize that there exists strong attraction in the range r < 1 
fm for the D*D* system with J = 0, 1 as can be seen clearly from the shape of the total effective potential in 
Fig.[5](a)-(b). 

If this enhancement is a molecular state, the potential candidate of Z+(4051) is the $**+ state with isospin 
/ = 1. Our numerical results indicate that there exists a bound state solution (1) for ^**+ with = 0"*" and 
A ~ 4 GcV; (2) for $**+ with — 1+ and A ~ 10 GeV. Unfortunately such a cutoff seems too large according 
to our criteria. 

If future experiments confirm Z"'"(4051) as a loosely bound molecular state, its quantum number is very probably 
= 0"*", which can be measured experimentally through the angular momentum distribution of tt'^XcI- Since 
^**+ belongs to an isospin triplet, its neutral partner state <i)**o may be searched in the tt'^XcI channel. 

2. X(3872) 

In our previous work [l^, we explored the D^D**^ molecular state by considering the tt and a meson exchange 
potentials only. Our numerical results showed that it is impossible to form a D'^D*'^ molecular state with the 
coupling constants determined by experiment and a reasonable cutoff around 1 GeV [3. Unfortunately we 
missed a minus sign in the sigma meson exchange potential, although its contribution is not large. Moreover we 
didn't take into account the charged modes. 

In this work, we consider both the neutral and charged modes and include the exchange force from the tt, 77, 
(7, p and u) mesons. The resulting total effective potential is attractive as shown in Fig. [TO] (f). The vector 
meson exchange especially the rho meson exchange provides the additional strong attraction. The ^g*' state 
corresponds the X(3872). 
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V-V 


of ill 


A (GeV) 


E (McV) 


rrms (fm) 




- 


- 


_ 


^ S 7 ^ S 7 ^ S 


2.80 


-1.58 


2.57 




2.90 


-15.38 


1.15 




1.90 


-16.33 


1.27 




2.00 


-26.25 


1.03 




- 


- 


- 


-til 

-i>:i 


0.72 


-13.66 


1.65 




0.73 


-8.65 


1.93 




0.75 


-22.85 


1.24 




0.77 


-4.66 


2.24 




0.75 


-28.48 


1.15 




0.78 


-35.20 


1.05 




1.21 


-4.47 


2.29 




1.23 


-19.31 


1.26 








A (GeV) 


E (McV) 


frms (fm) 


A*± o^t' (S*^ 










1.60 


-6.24 


1.74 




1.70 


-44.08 


0.74 


ii J, _l_ ^-1,(1 


0.87 


-5.65 


1.12 




0.90 


-12.67 


1.03 












1.15 


-5.70 


1.88 




1.20 


-30.52 


0.96 




0.82 


-34.13 


1.07 




0.84 


-3.57 


2.43 







TABLE XII; The possible bound state solutions for the P-V system after increasing all the coupling constants by a factor of 
two. 



Recall that the sigma meson exchange force is weak and the pion meson exchange alone does not bind the D*D 
into a molecular state with a physical pionic coupling constant and a reasonable cutoff [l^ . In fact a bound state 
appears only with A ~ 6 GeV there. With the additional strong attraction from the vector meson exchange, our 
present numerical analysis shows that a molecular state exists with A ~ 0.55 GeV. If we arbitrarily increase all 
the coupling constants by a factor of two, there exists a molecular state with a slightly more reasonable A ~ 0.75 
GeV. If we reduce all the coupling constants by a factor two, there appears a bound state with A ~ 0.45 GeV. 
Such a cutoff value is below 1 GeV and smaller than the exchanged vector meson mass. 

The pionic coupling constant of the D*D system was extracted from the decay width of the D* meson exper- 
imentally. Its value is known quite reliably. In contrast, the vector coupling constant was estimated using the 
vector dominance model in Ref. [28l |. which may be overestimated according to a recent light cone QCD sum 
rule analysis [30| . If we fix the pionic and scalar coupling constants but reduce the vector coupling by a factor 
of two only, there exists a very reasonable bound state solution with E = —1.75 GeV and A = 1.4GeV as shown 
in Table IXVlII Now the value of its root-mean-square radius reaches 3.53 fm! In other words, such a molecular 
state is bound very loosely. 

In short summary. A" (3872) may be accommodated as a molecular state dynamically, although drawing a very 
definite conclusion is very difficult especially when other available experimental information is taken into account. 
In our calculation, we assumed the SU(2) symmetry and ignored the mass difference between the charged and 
neutral mesons. Moreover we ignored the possible S-D mixing effect completely since we impose the S-wave 
condition in the derivation of the potential. These approximations together with the sensitivity of the numerical 
results to the cutoff parameter should be considered in the future investigation. 



3. A:(3764) 
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V-V 




A (GeV) 


E (MeV) 


rrms (fm) 


y¥>,*i 

<P (p CD 




(5*zt 












$*±,$*" 


0.36 


-11.10 


2.07 




0.37 


-7.30 


2.38 


^ si 


0.39 


-17.22 


1.74 




0.41 


-9.31 


2.11 


<I>*o 


0.43 


-12.33 


1.78 




0.45 


-5.42 


2.32 




0.44 


-18.01 


1.68 




0.46 


-8.66 


2.11 




0.38 


-13.56 


1.90 




0.40 


-6.17 


2.41 




B-B' 


State 


A (GeV) 


E (MeV) 


rrins (fm) 


n»± 0*0 n«o 




Q»± o*" O*" 




o*± 0*0 








1^*0 






3.30 


-1.39 


2.71 




3.40 


-11.03 


1.41 




0.46 


-19.73 


1.50 




0.48 


-8.10 


1.78 







TABLE XIII: The possible bound state solutions for the P-V system after reducing all the coupling constants by a factor of 
two. 
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state 


remark 


state 


remark 
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^ s •> ^ S 1 ^ S 
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? 
? 




® 
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TABLE XIV: Summary of the T) ~T> and B ~ B systems. The symbols and Q denote "forbidden" and "allowed" respectively. 



The broad structure X(3764) is close to the DD threshold and was observed in e+e~ — > hadrons process. Its 
quantum number is = 1^^. In the present case, we need consider the P-wave $g by adding the centrifugal 
potential £(1 -|-£)/(2/ir^) with £ — 1 into the potential in Eq. p9)) . Here /i is the reduced mass of the system. In 

Fig. [121 we show the variation of V{r)j?^^^i + 3/(2/xr^) with several typical values of A. Our numerical analysis 
indicates that there does not exist a P-wave <I>g from such a potential. In other words, the structure X(3764) is 
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TABLE XV: Summary of the V - V* and B* - B* systems. 
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TABLE XVI: Summary of the V-V and B - B* systems. 



A 


E (MeV) 


Trms (fm) 


1.6 


-12.86 


1.47 


1.5 


-6.07 


1.98 


1.42 


-2.41 


3.09 


1.41 


-2.07 


3.30 


1.40 


-1.75 


3.53 



TABLE XVII: The bound state solutions for the D* — D system if we fix the pionic and scalar coupling constants and reduce 
only the vector coupling constant in Section Hill by a factor of two. 

not a DD molecular state. 
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FIG. 12: The dependence of the potential + 3/(2/ir'^) on A. The solid, dotted and dashed lines correspond to 

A = 1, 3, 4 GeV respectively. 
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